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Abstract 

Background: Complex life histories require adaptation of a single organism for multiple ecological niches. 
Transitions between life stages, however, may expose individuals to an increased risk of mortality, as the process of 
metamorphosis typically includes developmental stages that function relatively poorly in both the pre- and post- 
metamorphic habitat. We studied predator-mediated selection on tadpoles of the wood frog, Rana sylvatica, to 
identify this hypothesized period of differential predation risk and estimate its ontogenetic onset. We reared 
tadpoles in replicated mesocosms in the presence of the larval odonate Anax Junius, a known tadpole predator. 

Results: The probability of tadpole survival increased with increasing age and size, but declined steeply at the 
point in development where hind limbs began to erupt from the body wall. Selection gradient analyses indicate 
that natural selection favored tadpoles with short, deep tail fins. Tadpoles resorb their tails as they progress toward 
metamorphosis, which may have led to the observed decrease in survivorship. Path models revealed that selection 
acted directly on tail morphology, rather than through its indirect influence on swimming performance. 

Conclusions: This is consistent with the hypothesis that tail morphology influences predation rates by reducing 
the probability a predator strikes the head or body. 



Background 

Many organisms exploit different environments over the 
course of their life cycle. Perhaps the most extreme 
example of this shift in resource use is that which 
accompanies metamorphosis in animals with complex 
life cycles [1]. Complex life cycles - hereafter referring 
to organisms with at least two discrete post-embryonic 
life-stages [2,3] - are ubiquitous in animals, being 
expressed in at least 80% of all species [4,5]. They may 
evolve for several reasons, such as trophic switching or 
specialized dispersal/breeding forms [6]. The tradeoffs 
that accompany shifts in niche occupancy will typically 
be accompanied by divergent selective regimes and 
alternative adaptations. In part, this accounts for the 
large differences in morphology, physiology, behavior, 
and other aspects of the phenotype observed among life 
stages. Although dramatically divergent morphologies 
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among different life stages allow individuals to exploit 
multiple kinds of resources throughout ontogeny, com- 
plex life cycles also involve functional trade-offs and 
thereby create a new problem: how to optimize the tran- 
sition between life stages [7,8]. 

The challenge of adapting to multiple adaptive peaks 
can be partially resolved by genetic and developmental 
decoupling among life stages [3]. Nonetheless, it is often 
the case that genetic, developmental, and functional cor- 
relations persist across life stages (e.g., [9-14]). More- 
over, even if there is complete adaptive decoupling of 
divergent life stages, the transitional period between life 
stages is still likely to be a performance trough that 
exposes individuals to increased risks. Indeed, the more 
differentiated the life stages, the more intense the risks 
are likely to be. Metaphorically, the transition from juve- 
nile to adult may be viewed as movement between alter- 
native peaks on an individual's adaptive landscape 
[15,16], where peaks represent correspondence between 
an individual's phenotypic traits and the local maximum 
probability of survival. 
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Many amphibians exhibit a complex life cycle in 
which larval development (intervals of which are 
referred to as Gosner stages in frog tadpoles; [17]) is fol- 
lowed by metamorphosis into an adult form [5,18]. Tad- 
poles are highly specialized for feeding, and the tadpole 
body plan consists mostly of a globose body and a 
sheet-like, laterally compressed tail [19]. During meta- 
morphosis, the tail is resorbed as hind and forelimbs 
emerge, thereby facilitating the transition from an aqua- 
tic swimming form (undulatory, axial locomotion) to a 
terrestrial hopping form (saltatory, appendicular loco- 
motion). It is during the intermediate stages of meta- 
morphosis that individuals are thought to experience 
increased predation risk [7]. The hypothesized period of 
increased predation risk separating larval and adult 
forms derives from the observation that metamorphs are 
optimized for neither larval nor adult niches [20]. For 
example, emergent hind limbs may impose drag and 
reduce swimming performance [19-21] and residual tail 
tissue may negatively impact saltatory locomotion [7]. 
For instance, metamorphosing chorus frogs, Pseudacris 
triseriata, are more likely to be captured by predatory 
garter snakes, Thamnophus sirtalis, than are either tad- 
poles or adult frogs [7]. Laboratory selection experi- 
ments on tadpoles likewise suggest the presence of a 
performance decline at metamorphosis [22]. Field 
experiments designed to measure both natural selection 
and variation in viability during ontogenetic stages near 
the developmental switch between life-stages, however, 
are still lacking. 

Here, we test a set of related hypotheses about varia- 
tion in survival probability in the wood frog, Rana sylva- 
tica. We begin broadly, by first testing whether fitness (i. 
e., survival) correlates with morphology across tadpole 
development [22-25]. This first question is designed to 
test the hypothesis of increased predation risk during 
tadpole metamorphosis. We next use path analytic mod- 
els to compare alternative hypotheses regarding the cau- 
sal structure underlying selection on tadpole 
morphology. This includes a test of the hypothesis that 
tail morphology is subject to selection via its effect on 
swimming performance [24,26], which may be important 
for predator escape. We also consider the alternative 
hypothesis that tail shape may, as has been demon- 
strated previously, enhance survival by serving as a 
"lure" to attract predatory attacks towards the tail, 
thereby reducing the probability of mortal wounds to 
the head/body region [25,27-30]. 

Methods 

We collected tadpoles of the wood frog, Rana sylvatica, 
from a single pond near Randolph, Vermont, USA (43° 
54' N, 72°38' W) on June 11, 2010. The pool naturally 
contained larval odonates and other predatory 



invertebrates (Calsbeek and Kuchta, pers. obs.). Tad- 
poles were held overnight in 5 gallon food-grade plastic 
buckets with filtered pond water, and were fed an ad 
libitum diet of boiled lettuce. The morning after cap- 
ture, we individually marked each tadpole with a unique 
color-coded combination of elastomer dyes (visible elas- 
tomer implants available from Northwest Marine Tech- 
nologies, Shaw Island WA, U.S.A.) that we injected into 
the dorsal half of the tail fin, posterior to the body wall. 
Tadpoles were immobilized (but not anesthetized) dur- 
ing the marking procedure by holding them in a plastic 
multi-channel pipette well. We scored each tadpole's 
developmental stage [17] with the aid of a dissecting 
microscope just prior to the initiation of the selection 
experiment (mean Gosner stage = 34 ± 4.62 SD). 

Tadpoles were then individually transferred to a V- 
shaped glass tank (which imposed a consistent orienta- 
tion on the tadpoles) with a size standard, and were 
digitally photographed. We used digital images of each 
tadpole to make the following linear measurements: 
Head length: the distance from the anterior tip of the 
snout to the junction of the body with the tail; head 
height: the depth of head at its tallest point; tail length: 
from the junction of the tail with the body wall to the 
distal tip of the tail; tail muscle height: muscle height at 
the tallest point of the tail muscle; and tail fin height: 
fin height at the tallest part of the tail. 

We measured swimming performance for half of the 
individuals in our selection study (N = 200 tadpoles) 
using a small (36L x 26W x 5H cm) tank containing fil- 
tered pond water and a size standard. Rapid development 
among the tadpoles held in buckets prevented us from 
measuring swimming performance for the remaining 200 
individuals. Swimming trials were videotaped at 250 
frames/sec using a high-definition digital camcorder (JVC 
Evario GZ-HM550-bu). Each tadpole was introduced to 
the swimming chamber and then motivated to initiate a 
"C-start" by touching the junction point between the tail 
fin and the body wall using a small metal pointer. We 
recorded three C-starts for each tadpole and used the 
fastest of these trials to estimate swimming performance, 
recording average speed over the 50 fastest frames. We 
chose to use this measure in our selection analyses 
because fifty frames was the average time required to 
swim one body length, and we assume that this is a good 
metric for predator avoidance. Burst speed was measured 
along the path of the tadpole movement using MotionA- 
nalysis software (available from M. Chappell, University 
of California, Riverside, CA, U.S.A. http://warthog.ucr. 
edu/). We used the tadpole eye as a landmark for track- 
ing individuals. All capturing, marking, photography, and 
swimming performance trials were conducted within 36 
hours and the tadpoles were immediately transferred to 
cattle tanks for the selection experiment. 
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Selection experiment 

We conducted our selection experiment using eight 
1136 L (300 gallon) cattle tanks that were randomly 
selected from an array of 49 tanks housed in an open 
field near the Dartmouth College campus. One month 
before introducing tadpoles, cattle tanks were cleaned 
and filled with ground water, 0.550 kg of dried Oak leaf 
litter, 15.4 g of rabbit chow (a nutrient source), and a 
three-liter aliquot of mixed zooplankton and phyto- 
plankton collected from a pond near Norwich, VT 
(43.73° N, 72.31°W). We added five larval dragonflies 
[Anax junius) to each tank to serve as predators on tad- 
poles. To provide developing frogs with retreat sites, we 
placed three White Water Lilly {Nymphea ordorata) 
fronds on the water surface of each tank. To provide 
dragonfly larvae with perches, we used stones to anchor 
three to five tree branches (~100 cm) to the bottom of 
each tank. Finally, we randomly assigned 50 tadpoles to 
each of the eight tanks. Set up this way, the cattle tanks 
functioned as self-sustaining mesocosms that mimicked 
conditions experienced by tadpoles in nature [31]. We 
covered each cattle tank with 0.5 x 0.5 cm hardware 
cloth pulled taut and secured with elastic cords. This 
functioned to shade the tanks, prevent predation, ensure 
that no metamorphosing individuals escaped, and pre- 
clude large, predatory insects from laying eggs. 

We recorded the identity of individual surviving tad- 
poles in our selection experiment five and fourteen days 
following introduction to the artificial ponds. Survival 
was scored by removing all the leaf litter and filtering 
each tank to recover tadpoles with hand-held dip nets. 
We also verified the presence of all five dragonfly-larvae 
in each mesocosm (one dragonfly-larva in each of three 
tanks was replaced to account for single dead indivi- 
duals). After the first census, we replaced the leaf litter, 
dragonfly-larvae, and tadpoles, and re-covered the tanks 
with the shade cloth. Following the second census, all 
tadpoles were brought back to the laboratory, sacrificed 
with an overdose of MS-222, and stored in 70% ethanol. 

General linear models were used to calculate selection 
gradients [32,33] for linear (P) and non-linear (y) forms 
of selection. First, competing models were compared 
using Akaike's information criterion (AIC, [34] Table 1). 
This metric, which does not require nested models, cal- 
culates the likelihood of the model given the data and 
the number of parameters. Consequently, two models 
with equal likelihoods that differ in the number of para- 
meters will have different AIC values, and the model 
with the smaller number of parameters will be favored. 
Next, the difference between the preferred model (i.e., 
with the lowest AIC value) and each of the subsequent 
models was calculated (At) by stepwise inclusion of the 
remaining linear and quadratic terms. We then calcu- 
lated the normalized relative likelihoods of the models, 



Table 1 Comparison of alternative selection models 
showing the number of parameters (k) in each model. 



Traits in the 
model 


k 


AIC 


Ai 


Likelihoods 


Wi 


Evidence 
Ratio 


Preferred model 


7 


893.91 


0 


1 


0.6066 




Head Height 


8 


895.7 


1.79 


0.4099 


0.2478 


2.4473 


Tail Muscle 
Height 


9 


897.69 


3.78 


0.151 


0.0916 


6.6194 


Head Height 2 


10 


899.81 


5.9 


0.0523 


0.0317 


19.106 


Tail Muscle 
Height 2 


11 


901.23 


7.32 


0.0257 


0.0156 


38.8613 


Tail Fin Height 2 


12 


902.95 


9.04 


0.0108 


0.0066 


91.8356 



AIC score = Akaike Information Criterion; A, = difference in AIC scores 
between the best model and subsequent model; w, = normalized relative 
likelihood of the model given the data; Evidence Ratio = the relative odds 
that a model is the best given the data. The first row of the table represents 
the preferred model, which includes linear terms for Head Length, Tail Length, 
and Tailfin Height, and quadratic terms for Head Length and Tail Length (see 
Table 4). Subsequent rows show the influence of adding each indicated trait 
in succession. Traits labeled with a superscript "2" represent quadratic terms. 

also known as the Akaike weights (wi), which quantify 
the relative support for different models [35]. Finally, we 
calculated the evidence ratio, which compares each 
model to the best model and provides the relative odds 
of competing models. 

Linear selection gradients were calculated from mod- 
els that included only linear terms, whereas quadratic 
gradients (e.g., non-linear selection) and cross-product 
terms (i.e., correlational selection) were calculated from 
models that included both the linear and quadratic 
terms. Because the GLM underestimates quadratic 
terms by half [36,37], quadratic gradients and their stan- 
dard errors were doubled. Though parametric statistics 
provide robust estimates of selection gradients and 
other parameters [32,38], these tests may be violated by 
survival data (live/die), which tend to have non-normally 
distributed errors [39,40]. We computed significance 
values for selection gradients using generalized linear 
models including a logit link function [41]. Prior to 
pooling data from individual tanks (i.e., replicates), we 
tested for any interaction between relevant terms and 
the factor for tank. None of these were significant, indi- 
cating that selection operated in the same way in all 
replicates. We dropped the interaction terms but 
retained a factor for "tank" in our models. The factor 
for tank explained a significant portion of the variance 
in all full models (0.02 > P < 0.03), but not in reduced 
models (0.06 > P < 0.08). All variables used in selection 
analyses were standardized to a mean of zero with unit 
standard deviation, except our fitness variable (survival), 
which was scaled by the mean [32,42]. The degree of 
multi-colinearity among traits was assessed by estimat- 
ing variance inflation factors (VIF; [43]), all of which 
were less than five. We visualized fitness surfaces using 
cubic splines [44]. 
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Path analysis was used to investigate the structure of 
causal relationships in our selection experiments [45,46]. 
First, we developed a set of a priori causal path models 
based on the competing hypotheses that variation in 
survival was most highly dependent on swimming per- 
formance versus predator evasion by caudal luring (Fig- 
ure 1). The swimming performance model included 
causal paths from morphology — > performance — > fit- 
ness, while the caudal luring model was reduced to cau- 
sal paths from morphology — > fitness only. This latter 
case models a situation in which a trait other than burst 
swim speed mediates the relationship between morphol- 
ogy and fitness [47]. A third model combined the two 
models above, and allowed for the possibility that mor- 
phology impacts fitness through both measured and 
unmeasured performance variables. Significance tests for 
individual path models were based on comparisons in 
which the covariance structure of each model was tested 
against the covariance expected under the assumption 
the model was correct [48]. A significant difference in 
this comparison indicates that the model in question 
provides a poor fit to the data. Path analyses, including 
significance tests, were performed using the program 
AMOS v. 18 [49] 

In a second analysis, our Model 3 (Figure 1) was itera- 
tively reduced to its significant components by sequen- 
tially setting causal paths with the lowest partial 
regression coefficients and the highest P values to zero 
[50]. Alternative models were compared using Akaike's 
information criteria (AIC, [34]), including the difference 
between the preferred model and each subsequent 
model (Ai), normalized relative likelihoods and evi- 
dence ratios [35,51]. 

In addition to path analysis using maximum likeli- 
hood, we also conducted Bayesian analyses of the data. 
We did this to account for the binomial distribution of 
our fitness variable (survival), which likely violates the 



assumption of normal errors and multivariate normality 
in least-square calculations [52,53]. Bayesian analysis in 
Amos 18 [49] employs a Markov Chain Monte Carlo 
(MCMC) algorithm for estimating posterior distribu- 
tions, and properly accounts for the binomial status of 
our fitness variable. Parameter estimates were estimated 
from 150,000 generations following a burn-in of 500 
generations. Convergence of the MCMC algorithm was 
assessed using the convergence statistic developed by 
Gelman et al. [54] and implemented in AMOS [49]. The 
significance of parameter estimates was assessed using 
95% Bayesian credibility intervals. The results of the 
Bayesian analysis of the path coefficients were very simi- 
lar to the maximum likelihood estimates (data not 
shown) and will not be presented. 

Results 

Despite the presence of floating refugia in each tank, our 
artificial ponds turned up one metamorphosed frog that 
clearly drowned after failing to find a terrestrial refuge. 
It is likely that some fraction of the mortality that we 
attributed to predation was from metamorphs that 
drowned. On the other hand, we also recovered six fully 
metamorphosed frogs that survived to the end of the 
experiment. As a conservative approach to analyses, we 
broke our data up into two different data sets. The first 
dataset is based on all individuals in the study, and is 
referred to as the "Full" dataset; it is described through- 
out the rest of this paper. The second dataset excluded 
all individuals whose Gosner stage was > 39 at the 
experiment's outset; this is referred to as the "reduced" 
data set. Based on rates of development measured in 
our study, this reduced dataset increases the chance that 
most individuals were in an aquatic stage throughout 
much of the course of the experiment. Morphological 
measurements could not be made for a few individuals 
and the size of these data sets varies slightly (see table 



Model 1 



Model 2 



Model 3 



Morphology 



Performance 



Fitness 



Morphology 



Fitness 



Morphology 



Performance 



Fitness 



Figure 1 Path diagram of the relationships between morphology, swimming speed, and fitness Model 1 is a classic path analysis 
diagram with links going from morphology to performance to fitness [62]. In Model 2, links go directly from morphology to fitness; performance 
is omitted. This would be the case if morphology impacted fitness through means other than swimming speed, for example by acting as a 
caudal lure. Model 3 is the "full model" and includes links between morphology and performance as well as between morphology and fitness, as 
would be the case if swimming speed as well as other factors mediated the relationship between morphology and fitness. 
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legends for details). Results from the two sets of ana- 
lyses were qualitatively nearly identical (Tables 2, 3 and 
4). 

In the full data set, the mean percentage survival (± 
SE) in each tank to the first census period (5 days) was 
0.77 ± 0.02 (range 0.66-0.90). By the second census (14 
days), mean survival had decreased to 0.57 ± 0.02 (range 
0.48-0.66). Qualitatively, selection results during the two 
time periods were nearly identical (data not shown), but 
to maximize our power to detect selection, and to sim- 
plify the presentation of results, we use viability esti- 
mates from the second census as our measure of fitness. 
Frequent bite marks on the tails of surviving tadpoles 
suggest that dragonfly larvae were a key source of mor- 
tality in our study populations. We also recovered two 
complete tails during our census, with elastomer tags 
still intact, from tadpoles that did not survive. We con- 
clude that mortality in the selection replicates was lar- 
gely due to predation by dragonfly larvae. 

Variation in survival was strongly linked to ontogeny 
(Gosner stage) and favored tadpoles at intermediate 
stages of development (quadratic effect of Gosner stage 
ANOVA F 

1,385 - 63.83, P < 0001) with a decline in sur- 
vival probability starting, on average, around Gosner 
stage 37 (Figure 2, center panel). We therefore included 

Table 2 Linear (p*) and quadratic (y) selection on all 
morphological traits measured in this study in a data set 
that included all individuals (Full, N = 172), as well as in 
a second dataset from which we excluded all individuals 
whose Gosner stage was > 39 at the time of release were 
excluded (Reduced, N = 151). 







Directional 


Quadratic 






selection 


selection 


Data 


Trait 


P ± SE 


y± SE 


set 








Full 


Head Length 


0.08 ±0.16 


- 0.06 ± 0.16 




Head Height 


0.10 ± 0.12 


0.14 ± 0.14 




Tail Length 


- 0.16 ± 0.13 


0.004 ± 0.16 




Tail muscle 


- 0.14 ± 0.15 


- 0.10 ± 0.14 




Height 








Tail Fin Height 


0.22 ± 0.09* 


- 0.06 ± 0.1 1 




Swimming speed 


- 0.01 ± 0.07 


- 0.06 ± 0.10 


Reduced 


Head Length 


0.05 ± 0.19 


- 0.10 ± 0.17 




Head Height 


0.09 ± 0.13 


0.14 ± 0.15 




Tail Length 


- 0.17 ± 0.14 


0.04 ±0.19 




Tail muscle 


- 0.08 ± 0.16 


- 0.03 ± 0.17 




Height 








Tail Fin Height 


0.21 ± 0.10* 


- 0.05 ± 0.12 




Swimming speed 


- 0.05 ± 0.08 


- 0.10 ± 0.11 



Asterisks indicate significant (* P < 0.05) selection, as determined by a 
generalized linear model with a logit link function and survival (0 or 1) as the 
response variable. A term for Gosner stage was included in these models, but 
because development stage is not a "trait" on which one can measure 
selection, is not shown here. A factor for tank was also included. See text for 
details. 



Table 3 Linear ((3) and quadratic (y) selection on all 
morphological traits except burst swimming speed. 







Directional 
selection 


Quadratic 
selection 


Data 
set 


Trait 


P ± SE 


y± SE 


Full 


Head Length 


- 0.01 ± 0.11 


- 0.18 ± 0.10 




Head Height 


0.10 ± 0.09 


0.04 ± 0.09 




Tail Length 


- 0.19 ± 0.07* 


- 0.16 ± 0.06* 




Tail muscle 
Height 


0.02 ± 0.09 


- 0.03 ± 0.09 




Tail Fin Height 


0.24 ± 0.06*** 


- 0.07 ± 0.06 


Reduced 


Head Length 


- 0.08 ± 0.12 


- 0.24 ±0.11* 




Head Height 


0.07 ±0.10 


0.07 ± 0.10 




Tail Length 


- 0.14 ± 0.08 


- 0.1 1 ± 0.08 




Tail muscle 
Height 


0.07 ±0.10 


-0.03 ± 0.10 




Tail Fin Height 


0.16 ± 0.07* 


-0.08 ± 0.08 



The results of two analyses are presented: a Full dataset that included all 
individuals (N = 381), as well as in a second dataset in data set (Reduced, N = 
324), from which we excluded all individuals whose Gosner stage was > 39 at 
the time of release. Asterisks indicate significant (* P < 0.05, ***P < 0.005) 
selection, as determined by a generalized linear model with a logit link 
function and survival (0 or 1) as the response variable. A term for Gosner 
stage was included in these models, but because development stage is not a 
"trait" on which one can measure selection, is not shown here. A factor for 
tank was also included. See text for details. 

a term for Gosner stage in selection models. For com- 
pleteness, we present models that include all measured 
traits (Table 2), models without swimming speed (which 
maximizes our sample size; Table 3), and a model using 
the set of independent variables corresponding to the 
smallest AIC score (Table 4). This last model, which we 
consider the preferred model (Table 1), included a linear 
term for tail fin height and Gosner stage, and linear and 
quadratic terms for tail length and head length. In this 

Table 4 Linear ((3) and quadratic (y) selection on 
morphological traits used in a model chosen based on 
AIC scores. 



Directional selection Quadratic selection 



Data set 


Trait 


P± SE 


y± SE 


Full 


Head Length 


0.05 + 0.09 


- 0.20 ± 0.07* 




Tail Length 


- 0.17 ± 0.07* 


- 0.16 ± 0.06* 




Tail Fin Height 


0.27 ± 0.05*** 




Reduced 


Head Length 


- 0.007 ± 0.10 


- 0.20 ± 0.08* 




Tail Length 


- 0.1 1 ± 0.08 


- 0.10 ± 0.08 




Tail Fin Height 


0.20 ± 0.07*** 





Two identical analyses are presented: the first included all individuals (Full, N 
= 381), the second is from a censured dataset (Reduced, N = 325) from which 
all individuals whose Gosner stage was > 39 at the time of release were 
excluded. Asterisks indicate significant (* P < 0.05, ***P < 0.005) selection, as 
determined by a generalized linear model with a logit link function and 
survival (0 or 1) as the response variable. A term for Gosner stage was 
included in these models, but because development stage is not a "trait" on 
which one can measure selection, is not shown here. A factor for tank was 
also included. See text for details. 
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25 30 35 40 45 



25 30 35 40 45 



25 30 35 40 45 



Developmental stage (Gosner) 



Figure 2 The probability of survival was maximized for tadpoles in intermediate stages (e.g., Gosner stages 30-37) of development 

Variation in survival is shown for the eight selection replicates used in our experiment (outer panels) and for all replicates pooled (central panel). 
The dark line represents the best fit cubic spline and light lines indicate the 95% confidence limits. The central panel includes illustrations of 
representative tadpoles/metamorphs at several developmental stages (indicated as numeric Gosner stages above each illustration). 



model, selection favored individuals with deep tail fins 
(P = 0.27 ± 0.05, P < 0.0001) and short tails (p = -0.17 
± 0.07, P = 0.02) (Figure 3). We also detected quadratic 
components to selection on tail length and head length 
that were both stabilizing (tail length: y l:l , = -0.16 ± 
0.06, P = 0.01; head length: y 2 , 2 = -0.20 ± 0.07, P = 0.01) 
(Figure 3). To verify that the results were not biased by 
the relationship between size and Gosner stage, we 
regressed Gosner stage against tail and head morphol- 
ogy, and saved the residuals. Patterns of selection based 
on residual trait values were qualitatively similar to 
those using raw values (e.g., selection for deep residual 



tail fins [p = 0.22 ± 0.04, P < 0.0001] and short residual 
tail lengths [p = -0.10 ± 0.04, P = 0.03]). 

Swimming speed was positively correlated with devel- 
opmental stage (r 2 = 0.14, df = 186, P < 0.0001) and tail 
length (r 2 = 0.12; df = 186, P < 0.0001), and there was a 
weak quadratic relationship between swimming speed 
and tail fin height (individuals of intermediate tail fin 
height swam fastest: r 2 = 0.07, df = 186, P = 0.051). 
However, we did not detect any selection on swimming 
speed in our experiment. In a model that included linear 
terms for tail length and tail fin height, the selection 
gradient for swimming speed was weakly negative and 
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Figure 3 In the preferred model (Table 3) chosen based on AIC scores (Table 1), natural selection acted primarily on tail length, tail 
height and head. Selection on tail length and head length both were stabilizing around shorter values. Selection on tail height was directional 
and positive. See the text and Table 3 for statistical details. The dark line represents the best fit cubic spline for each trait, and light lines indicate 
the 95% confidence limits. 



non-significant (P = -0.03 ± 0.07, P = 0.65). Even when 
we considered selection on swimming speed alone (i.e., 
the selection differential for swimming speed) we 
detected no variation in survival that was related to 
swimming performance (s = -0.07 ± 0.06, P = 0.23). 

The results of the path analyses parallel the multiple 
regression analyses. The best fit model was Model 2 
(Morphology -> Fitness; AIC = 52.41; DIC = 1056.30; 
Figure 1; Table 5). Over 99% of the relative likelihood 
was captured by this model, and the relative odds of the 
second best model being better than the most strongly 
supported model was 3159:1 (Table 5). In contrast, 
Model 1 (Morphology — > Performance — > Fitness; Figure 
1) was significantly different from the data (% 2 = 46.88; 
P < 0.001). We thus conclude that Model 2 is strongly 
supported relative to alternatives. 

In our second path analytic approach, we iteratively 
reduced causal paths by removing the most poorly sup- 
ported paths after each run until we were left with a 
model in which all causal paths were significant. We 
started with the Full Model (Model 3; Figure 1) because 
this model included all theoretically interesting causal 
paths. The most fully reduced model received the stron- 
gest support (AIC = 61.38; Table 6). In this best fit 
model, the only significant causal paths were between 
tail fin height and fitness (P = 0.24; P < 0.001) and tail 
length and fitness (P = -0.32; P < 0.001). However, the 



best-fit model was not a robust improvement over 
related models. For example, the difference in AIC 
between the best-fit model and the 4 th best model was 
only 1.00, and for the 5 th best model, 2.19 (Table 6). In 
addition, the relative odds of models 2-4 ranged from 
1.35-1.65:1, and the relative odds of model 5 was 3:1 
(Table 6). We conclude that the first five models are not 
easily distinguished, we therefore show the results of 
Model 5 in our path analysis diagram because it is the 
fullest model receiving statistical support (Figure 4). 
Relative to the best-fit mo del, the 5 th best model 
includes causal paths between head height — > fitness, 
head length — > fitness, tail muscle height — > maximum 
swim speed, and maximum swim speed — > fitness (Fig- 
ure 4). None of these were significant, however, there is 
a trend between head length — > fitness (P = -0.16; P = 
0.06). Note that the link between swimming perfor- 
mance — > fitness is weak and not significant (P = 0.04; P 
> 0.05). 

Discussion 

One of the most common hypotheses regarding the evo- 
lution of complex life cycles is that alternative morpho- 
logical strategies are employed to exploit different 
resources throughout ontogeny. The transition between 
life stages, however, can be a vulnerable period in which 
individuals suffer higher rates of mortality. We have 



Table 5 Comparison of path models. 

Maximum Likelihood Analyses 
Model % 2 df P AIC fl, w. Evidence Ratio 

Morphology - > Fitness 0.41 1 0.52 52.41 0 0.9997 

Full Model 0.53 1 0.47 68.53 16.12 3.16 X Iff 4 3158.97 

Morphology -> Performance -> Fitness 46.88 5 < 0.001 106.88 54.47 1.48 xlO" 12 6.74 x 10 11 

Chi-square values (x 2 ), degrees of freedom {df), and the associated P value report the significance of the model. AIC = Akaike Information Criterion; A, - 
difference in AIC scores between the best model and subsequent model; w t = normalized relative likelihood of the model given the data; Evidence Ratio = the 
relative odds that a model is the best given the data. 
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Table 6 Comparison of path models. 



Model 


Path constrained to zero 


x 2 


df 


P 


Maximum Likelihood Analyses 
AIC A, w. 


Evidence Ratio 


1 


Head length — > Fitness 


5.38 


1 


0.61 


61.38 


0 


0.28 




2 


Max. swim speed -> Fitness 


1.98 


5 


0.85 


61.98 


0.60 


0.21 


1.35 


3 


Head ht. -» Fitness 


4.28 


6 


0.64 


62.28 


0.89 


0.18 


1.56 


■'1 


Tail muscle ht. — > Max. swim speed 


8.39 


8 


0.40 


62.39 


1.00 


0.17 


1.65 


5 


Tail fin ht. — > Max. swim speed 


1.58 


A 


0.81 


63.58 


2.19 


0.09 


3.00 


6 


Tail length -» Max. swim speed 


1.20 


3 


0.75 


65.20 


3.81 


0.04 


6.73 


7 


Tail muscle ht. -» Fitness 


0.89 


2 


0.64 


66.89 


5.51 


0.02 


15.70 


8 


Head ht. — > Max. swim speed 


0.36 


1 


0.55 


68.36 


6.98 


0.01 


32.77 


9 


Saturated model 




0 




70.00 


8.62 


0.00 


74.33 



Chi-square values (x 2 ), degrees of freedom {df), and the associated P value report the significance of the model. Note that model 9 is saturated, and thus the fit 
of the model to the data could not be tested using the chi-square statistic. In models that are not significant, the data are a good fit to the model. AIC = Akaike 
Information Criterion; A/ = difference in AIC scores between the best model and subsequent model; w, - normalized relative likelihood of the model given the 
data; Evidence Ratio = the relative odds that a model is the best given the data. 



presented empirical evidence that tadpoles of the wood 
frog, Rana sylvatica, when facing predation by dragonfly 
larvae, experience a higher probability of mortality as 
they approach metamorphic climax. That mortality 



probabilities increase during metamorphosis is not 
unexpected, as a tadpole with emergent hind and fore- 
limbs is well adapted for neither swimming nor jumping 
[7,20,21]. For instance, Arnold and Wassersug [8] 




Figure 4 Path diagram of the relationships between morphology, swimming speed, and fitness. The results of model 5 (Table 5) are 
illustrated here. Double4neaded arrows represent covariances (range: 0.551 - 0.854), and all of them are significant. Values near single4ieaded 
arrows are maximum likelihood parameter estimates of partial regression coefficients (direct effects. Arrows lacking a number represent causal 
paths set to zero (Table 5). Arrow thickness is proportional to the strength of relationship. Black arrows present significant parameter estimates, 
and grey arrows represent relationships that are not significant. The P-value of the dark grey arrow (head length -» fitness) is 0.06. Note that the 
causal paths between head length — > fitness and head length -» maximum swim speed are significant in model 9, which had the best AIC 
score (Fable 6). 
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showed across a large geographic range (Mexico to 
Washington state) that garter snakes, Thamnophis spp., 
were more likely to have consumed anuran metamorphs 
(tree frogs and toads) than either tadpoles or adults. 
They concluded that transforming anurans were highly 
susceptible to snake predation as a consequence of 
"locomotor ineptitude." 

Our data further suggest that selection acts strongly 
on morphological traits, favoring tadpoles with short 
tails and deep tail fins, but that this selection acts largely 
independently of swimming performance. This latter 
result is surprising given that tail shape influences swim- 
ming performance [55]. Indeed, in our data swimming 
performance was correlated with both tail length and 
tail fin depth, and larger values of both tail elements 
produced greater swimming speeds, consistent with pat- 
terns demonstrated elsewhere [56,57]. Our analyses may 
have suffered from reduced power given that we could 
only measure swimming speed for half of our study ani- 
mals. However, even when we removed all other terms 
from the model and measured selection differentials on 
swimming speed alone, the results were not significant. 
Moreover, path analyses revealed that the effects of 
morphology (tail length and tail fin height) were largely 
direct, acting to enhance survival probability per se, 
rather than serving as a functional link to swimming 
performance. We interpret this result as consistent with 
the hypothesis that short tails and deep tail fins are 
adaptive because they attract predatory strikes and 
increase the probability that a predator will strike tail 
tissue rather than sites on the head or body (i.e., "the 
caudal lure hypothesis"; [25,27,30]). 

Tadpoles of many frog species exhibit developmental 
plasticity in response to chemical cues from potential 
predators, whereby they develop a relatively deep tail fin 
and a small body {e.g., [22,58-61]). In particular, 
enlarged tail fins lead to enhanced survival in the pre- 
sence of larval odonates (summarized in [58]). There is 
reason to believe, however, that differences in tail shape 
do not influence swimming performance effectively 
enough to have a large impact on survival in the pre- 
sence of odonate larvae. This result is unexpected at 
first blush, given the high prevalence of causal relation- 
ships between morphology and performance in other 
animal systems [62-64]. Van Buskirk and McCollum 
[24] used experimental manipulation of tail fin morphol- 
ogy, trimming tissue to reduce both the total length and 
depth of the tail fin, to investigate the direct effects of 
changes in tail morphology on swimming performance. 
Their study revealed that changes in swimming perfor- 
mance were not apparent until one third of the tail was 
surgically removed, leading them to conclude that 
reduced susceptibility to predation must have been due 
to something other than enhanced swimming 



performance. Similarly, Wilbur and Semlitch [65] 
showed that damaged tails of R. utricularia incurred lit- 
tle survivorship cost in the presence of predatory newts 
{Notophthalmus viridescens). On the other hand, Van 
Buskirk et al [28] showed that tadpoles with predator- 
induced morphologies suffered fewer lethal strikes to 
the body, suggesting that enlarged tail fins may enhance 
survival via a "caudal lure" effect. 

The approach adopted in this study was to quantify 
relative survival and selection across ontogeny. One 
challenge faced by such an approach is that changes in 
size and shape are confounded throughout the develop- 
ment of the tadpole. This is the phenomenon summar- 
ized by Gosner stages. In addition, we were only able to 
quantify swimming performance and morphometric 
variables at the start of the study. Depredated tadpoles, 
unfortunately, cannot be measured. Our analyses thus 
assume that fundamental elements of size and shape 
were captured in our initial measures, and that the 
quantitative signal is maintained to some degree 
throughout ontogeny. If this were not the case, it is 
unlikely that we would have obtained sensible results. 

Though the number of studies of natural selection has 
grown rapidly in recent decades [66,67], there are still 
fundamental gaps in our understanding of the selective 
process. This is, in part, owing to the fact that selection 
studies are rarely replicated either temporally or spatially 
[41] and when studies are replicated, selection estimates 
tend to be highly variable among replicates [68]. Our 
study provides a rare example of repeatable selection, as 
replicate estimates of selection were highly congruent 
among mesocosms, suggesting that the changes that 
characterize metamorphosis are subject to strong and 
consistent patterns of selection among individuals. 

Conclusions 

Our study demonstrates an increase in mortality risk as 
tadpoles began to metamorphose. Owing to the nature 
of our experimental design, which focused on tadpole 
mortality, our data did not examine the effects of the 
transition from tadpole to froglet on survivorship (see 
[8]). As metamorphosis proceeds and the tail fin is 
resorbed, we expect that froglets would become better 
at hopping and thus less susceptible to predation. We 
suggest, as have others [19], that selection should thus 
favor individuals that minimize the transition time dur- 
ing metamorphic climax. This does not necessarily 
mean that selection should favor the most rapid possible 
development. Indeed, faster overall development often 
results in small adult body sizes, a condition that can 
have serious fitness consequences for adult anurans 
[69,70]. Rather, the optimal strategy should be to meta- 
morphose at a rate that maximizes the balance between 
the probability of surviving metamorphosis and later 
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fitness costs. Future studies should aim to measure 
selection on the separate components of developmental 
timing to improve our understanding of the targets of 
selection, including the costs and benefits of pursuing 
alternative metamorphic strategies. 

Acknowledgements 

We thank David Herron, Ariel le Filiberti, and Kelsey Pinson for help in the 
field. Thanks to Robert Cox and two anonymous referees for helpful 
suggestions that improved the manuscript. All research was conducted 
under iACUC protocol # 10-02-01 and funding was provided by Dartmouth 
College. RC was also supported by the National Evolutionary Synthesis 
Center. 

Author details 

'Department of Biological Sciences, Dartmouth College, Hanover, NH 03755 
USA. 2 National Evolutionary Synthesis Center, Durham, NC 27606 USA. 
department of Biological Sciences, Ohio University, Athens, OH 45701 USA. 

Authors' contributions 

RC and SK designed and carried out the experiments, analyzed the data, 
and wrote the manuscript. Both authors read and approved the final 
manuscript. 

Received: 10 August 2011 Accepted: 7 December 2011 
Published: 7 December 2011 

References 

1. Ebenman B: Evolution in organisms that change their niches during the 
life-cycle. American Naturalist 1992, 139(5)390-1021. 

2. Wilbur HM: Complex life-cycles. Annual Review of Ecology and Systematics 
1980, 11:67-93. 

3. Moran NA: Adaptation and constraint in the complex life-cycles of 
animals. Annual Review of Ecology and Systematics 1994, 25:573-600. 

4. Werner EE, Hall DJ: Ontogenetic habitat shifts in bluegill: The foraging 
rate-predation risk trade-off. Ecology 1988, 69:1352-1366. 

5. Hanken J, Hall BK, Wake MH: Larvae in amphibian development and 
evolution. The origin and evolution of larval forms 1999, 61-108. 

6. Istock CA: Evolution of complex life cycle phenomena- An ecological 
perspective. Evolution 1967, 21(3)592-605. 

7. Wassersug RJ, Sperry DG: Relationship of locomotion to differntial 
predatin on Pseudacric-tridseriata (Anura Hylidae). Ecology 1977, 
58(4):830-839. 

8. Arnold SJ, Wassersug RJ: Differential predation on metamorphic anurans 
by gartner sankes (Thamnophis)- Social behavior as a possible defense. 

Ecology 1978, 59(5):1 014-1 022. 

9. Shaffer HB, Austin CC, Huey RB: The consequences of metamorphosis on 
salamander (Ambystoma) locomotor performance. Physiological Zoology 

1991, 64(1)212-231. 

10. Blouin MS: Genetic correlations among morphometric traits and rates of 
growth and differentiation in the green tree frog, Hyla cinerea. Evolution 

1992, 46(3):735-744. 

11. Van Buskirk J, Saxer G: Delayed costs of an induced defense in tadpoles? 
Morphology, hopping, and development rate at metamorphosis. 
Evolution 2001, 55(4):821-829. 

12. Watkins TB: A quantitative genetic test of adaptive decoupling across 
metamorphosis for locomotor and life-history traits in the Pacific tree 
frog, Hyla regilla. Evolution 2001, 55(8):1 668-1 677. 

13. Azizi E, Landberg T: Effects of metamorphosis on the aquatic escape 
response of the two-lined salamander (Eurycea bislineata). Journal of 
Experimental Biology 2002, 205(6)841-849. 

14. Wilson RS, Kraft PG, Van Damme R: Predator-specific changes in the 
morphology and swimming performance of larval Rana lessonae. 
Functional Ecology 2005, 19(2):238-244. 

15. Dejongh HJ: Functional morphology of the jaw aparatus of larval and 
metamorphosing Rana temporaria. Netherlands Journal of Zoology 1968, 
18(1)1-103. 



16. Svensson El, Calsbeek R, (eds.): Evolution on Adaptive Landscapes. Oxford, 
UK: Oxford University Press: 2012. 

17. Gosner KL: A simplified table for staging anuran embryos and larvae 
with notes on identification. Herpetologiea 1960, 16:183-190. 

18. McDiarmid RW, Altig R: Tadpoles: The biology of anuran larvae. Tadpoles: 
The biology of anuran larvae 1 999, i-xiv, 1 -444. 

19. Wassersug RJ: Locomotion in amphibian larvae (or why aren't tadpoles 
built like fishes). American Zoologist 1989, 29(1)65-84. 

20. Huey RB: Sprint velocity of tadpoles (Bufo boreas) through 
metamorphosis. Copeia 1980, , 3: 537-540. 

21. Dudley R, King VA, Wassersug RJ: The implications of shape and 
metamorphosis for drag forces on a generalized pond tadpole (Rana 
catesbeiana). Copeia 1991, , 1: 252-257. 

22. Van Buskirk J, Relyea RA: Selection for phenotypic plasticity in Rana 
sylvatica tadpoles. Biol J Linn Soc 1998, 65(3)301-328. 

23. Van Buskirk J, McCollum SA, Werner EE: Natural selection for 
environmentally induced phenotypes in tadpoles. Evolution 1997, 
51:1983-1992. 

24. Van Buskirk J, McCollum SA: Influence of tail shape on tadpole swimming 
performance. Journal of Experimental Biology 2000, 203(14)2149-2158. 

25. Johnson JB, Burt DB, DeWitt TJ: Form, function, and fitness: Pathways to 
survival. Evolution 2008, 62(5):1 243-1 251. 

26. Teplitsky c, Plenet S, Lena JP, Mermet N, Malet E, Joly P: Escape behaviour 
and ultimate causes of specific induced defences in an anuran tadpole. 
Journal of Evolutionary Biology 2005, 18:180-190. 

27. Caldwell JP: Disruptive selection- a tail color polymorphoism in Acris 
tadpoles in response to differential predation. Canadian Journal of 
Zoology-Revue Canadienne De Zoologie 1982, 60(1 1)2818-2827. 

28. Van Buskirk J, Anderwald P, Lupoid S, Reinhardt L, Schuler H: The lure 
effect, tadpole tail shape, and the target of dragonfly strikes. Journal of 
Herpetology 2003, 37(2):420-424. 

29. Van Buskirk J, Aschwanden J, Buckelmuller I, Reolon S, Ruttiman S: Bold tail 
coloration protects tadpoles from dragonfly strikes. Copeia 2004, , 3: 
599-602. 

30. Johnson KE, Eidietis L: Tadpole body zones differ with regard to strike 
frequencies and kill rates by dragonfly naiads. Copeia 2005, , 4: 909-913. 

31. Semlitsch RD, Boone MD, Dodd CK Jr: Aquatic mesocosms. Amphibian 
ecology and conservation: a handbook of technigues [Technigues in Ecology & 
Conservation Series] 2010, 87-104. 

32. Lande R, Arnold SJ: The measurement of selection on correlated 
characters. Evolution 1983, 37:1210-1226. 

33. Brodie EDI, Moore AJ, Janzen FJ: Visualizing and quantifying natural 
selection. Trends Ecol Evol 1995, 10(8)313-318. 

34. Akaike H: Factor-Analysis and Aic. Psychometrika 1987, 52(3)317-332. 

35. Burnham KP, Anderson DR: Model selection and multimodel inference: a 
practical information-theoretical approach. 2002. 

36. Phillips PC, Arnold SJ: Visualizing Multivariate Selection. Evolution 1989, 
43(6):1 209-1 222. 

37. Stinchcombe JR, Agrawal AF, Hohenlohe PA, Arnold SJ, Blows MW: 
Estimating nonlinear selection gradients using quadratic regression 
coefficients: Double or nothing? Evolution 2008, 62(9)2435-2440. 

38. Mitchell-Olds T, Shaw RG: Regression analysis of natural selection: 
statistical and biological interpretation. Evolution 1987,41:1149-1161. 

39. Brodie ED III, Janzen FJ: On the assignment of fitness values in statistical 
analyses of selection. Evolution 1996, 50:437-442. 

40. Janzen FJ, Stern HS: Logistic regression for empirical studies of 
multivariate selection. Evolution 1998, 52(6): 1 564-1 571. 

41 . Svensson El, Eroukhmanoff F, Friberg M: Effects of natural and sexual 
selection on adaptive population divergence and premating isolation in 
a damselfly. Evolution 2006, 60(6):1 242-1 253. 

42. Brodie ED: Correlational Selection for Color Pattern and Antipredator 
Behavior in the Garter Snake Thamnophis ordinoides. Evolution 1992, 
46(5):1 284-1 298. 

43. Petraitis PS, Dunham AE, Niewiarowski PH: Inferring multiple causality: The 
limitations of path analysis. Functional Ecology 1996, 10(421-431). 

44. Schluter D, Nychka D: Exploring fitness surfaces. Amer Nat 1994, 
143(4)597-616. 

45. Wright S: The method of path coefficients. Ann Math Stats 1934, 
5:161-215. 

46. Li CC: Path Analysis - a primer. Pacific Grove, California: Boxwood Press; 
1975. 



Calsbeek and Kuchta BMC Evolutionary Biology 201 1, 1 1:353 
http://www.biomedcentral.eom/1 471 -2 1 48/1 1 /353 



Page 11 of 1 1 



47. Garland T Jr, Losos JB: Ecological morphology of locomotor performance 
in squamate reptiles. In Ecological Morphology: Integrative Organismal 
Biology Edited by: Wainwright PC, Reilly SM. Chicago: University of Chicago 
Press: 1994:. 

48. Gomez JM, Zamora R: Spatial variation in the selective scenarios of 
Hormatophylla spinosa (Cruciferae). American Naturalist 2000, 155:630-648. 

49. Arbuckle JL: Amos 18 User's Guide. Chicago, IL: SPSS Inc.: 2009. 

50. Gomez JM, Perfectti F, Camacho JPM: Natural selection on Erysimum 
mediohispanicum flower shape: insights into the evolution of 
zygomorphy. American Naturalist 2006, 168:531-545. 

51. Carstens BC, Stoute HN, Reid NM: An information-theoretic approach to 
phylogeography. Molecular Ecology 2009, 18:4270-4282. 

52. Tabachnick BG, Fidell LS: Using multivariate statistics. New York, NY: 
HarperCollins College:, 3 1996. 

53. Arbuckle JL: AMOS - Analysis of Moment Structures. Psychometrika 1994, 
59(1):1 35-137. 

54. Gelman A, Carlin JB, Stern HS, Rubin DB: Bayesian Data Analysis. Boca 
Raton: Chapman and Hall/CRQ, 2 2004. 

55. Dayton GH, Saenz D, Baum KA, Langerhans RB, DeWitt TJ: Body shape, 
burst speed and escape behavior of larval anurans. Oikos 2005, 
111(3)582-591. 

56. Wassersug RJ, Hoff K: The kinematics of swimming in anuran larvae. 
Journal of Experimental Biology 1985, 119:1-30. 

57. Liu H, Wassersug RJ, Kawachi K: A computational fluid dynamics study of 
tadpole swimming. Journal of Experimental Biology 1996, 199(6):1 245-1 260. 

58. Van Buskirk J: Natural variation in morphology of larval amphibians: 
Phenotypic plasticity in nature? Ecological Monographs 2009, 
79(4):681-705. 

59. Van Buskirk J, McCollum SA: Functional mechanisms of an inducible 
defence in tadpoles: morphology and behaviour influence mortality risk 
from predation. Journal of Evolutionary Biology 2000, 13(2)336-347. 

60. Relyea RA: The many faces of predation: How induction, selection, and 
thinning combine to alter prey phenotypes. Ecology 2002, 

83(7):1 953-1 964. 

61. Relyea RA: Fine-tuned phenotypes: Tadpole plasticity under 16 
combinations of predators and competitors. Ecology 2004, 85(1 ):1 72-1 79. 

62. Arnold SJ: Morphology, performance and fitness. Am Zool 1983, 
23:347-361 . 

63. Wainwright PC, Reilly SM: Ecological Morphology: Integrative Organismal 
Biology. Chicago, ILL: University of Chicago Press: 1999. 

64. Schwenk K: Feeding: Form, Function, and Evolution in Tetrapod 
Vertebrates. San Diego, CA: Academic Press: 2000. 

65. Wilbur HM, Semlitsch RD: Ecological consequences of tail injury in Rana 
tadpoles. Copeia 1990, , 1: 18-24. 

66. Kingsolver JGH, Hoekstra HE, Hoekstra JM, Berrigan D, Vignieri SN, Hill CE, 
Hoang A, Gilbert P, Beerli P: The strength of phenotypic selection in 
natural populations. Am Nat 2001, 157:245-261. 

67. Hereford J, Hansen TF, Houle D: Comparing strengths of directional 
selection: how strong is strong? Evolution 2004, 58:2133-2143. 

68. Siepielski AM, DiBattista JD, Carlson SM: It's about time: the temporal 
dynamics of phenotypic selection in the wild. Ecology Letters 2009, 
12(1 1):1 261 -1276. 

69. Merila J, Laurila A, Lindgren B: Variation in the degree and costs of 
adaptive phenotypic plasticity among Rana temporaria populations. 

Journal of Evolutionary Biology 2004, 17(5): 1 132-1 140. 

70. Altwegg R, Reyer HU: Patterns of natural selection on size at 
metamorphosis in water frogs. Evolution 2003, 57(4):872-882. 



doi:10.1 186/1471-2148-1 1-353 

Cite this article as: Calsbeek and Kuchta: Predator mediated selection 
and the impact of developmental stage on viability in wood frog 
tadpoles (Rana sylvatica). BMC Evolutionary Biology 201 1 11 :353. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f\ RioM _-| rpntr ,i 

www.biomedcentral.com/submit \ J ™°l™a central 



